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Abstract. We construct some intrinsically denned discrete model of the mag- 
netic Laplacian. The existence and uniqueness of solutions of the Dirichlct 
problem for the difference Poisson type equation are proved. We study in 
detail properties of the discrete model including the limiting process in the 
two-dimensional Euclidean case. 

1. Introduction 

Let (M,g) be a Riemannian manifold with a Riemannian metric (<?y) and 
dimM = n. Denote A P {M) the set of all diffcrcntiablc complex- valued p- forms on 
M for each p = 0, 1, ...,n. Note that A°(M) is just C°°(M). Denote also A? fe) (M) 

the set of all fc-smooth (of the class C k ) complex- valued p- forms on M . We define 
a magnetic potential as a real- valued 1-form A G A^(M). So in local coordinates 
x\,...,x n it can be written as 

n 

A = ^Ajdx\ 

3 = 1 

where Aj = Aj(x) are real- valued functions of the class C 1 . 

Let * be the metric adjoint operator (Hodge star operator) * : A P (M) — * 
A n ~ p (M). Then the invariant inner product of p- forms with compact support 
can be defined in a standard way by the relation 

(1.1) (<p, ip) = J (pA*ij, 

M 

where the bar over ip means the complex conjugation and A is the exterior multipli- 
cation operation. It is known that using the inner product in spaces of smooth 
p-forms with compact support we can define the completions of these spaces. We 
denote these Hilbert spaces by L 2 (M) for 0-forms (functions) and by L 2 A P (M) for 
p- forms, p = 1, n. Let us define the operator 

d: L 2 A p -\M) -> L 2 A p (M) 

as the closure in the i 2 -norm generated by the inner product HI. If) of the corre- 
sponding operator specified on smooth forms, i.e. as a strong extensions of the 
differential operator d : A P ~ 1 (M) — > A P (M). We will need a deformed differential 

(1.2) d A : L 2 {M) -> I?A X {M), tp^dtp + itpA, 
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where i = y— T and A is the magnetic potential. 

The definition of the invariant inner product immediately induces the formal 
adjoint operator to the differential operator cLa- So we have the operator 

5 A ■ L 2 A 1 (M) L 2 (M) 

giving by the relation 

(d A <p, u) = (<p, S A u), if G L 2 (M), lo G L 2 A 1 (M). 

Here we assume that one of the forms have compact support. Then we can define 
the magnetic Laplacian as follows 

(1.3) -A A = S A d A : L 2 (M) -> L 2 (M). 

Let us identify the magnetic potential A with the multiplication operation 

(1.4) A : L 2 (M) — ► i 2 A 1 (M), tp-upA 

(see e.g. |Jj). Then the formally adjoint operator 5^ can be written as follows 

(1.5) 8 A uj = {8-iA*)uj, 

where 5, A* are the formal adjoint operators to d and A respectively. Using 
l|1.2[l . l|1.3[l . we can rewrite the magnetic Laplacian as follows 

-A A tp = (S- iA*){d(f + iA<p) 

= 6dip - iA*dip + iS(Aip) + A*Aip 
= -Aip - iA*dip + i6(Acp) + A*Acp. 

Operator (|1.3|) is essentially self-adjoint (see for details J7| Th. 6.1]). 

The main purpose of this paper is to construct an intrinsically defined discrete 
model of the magnetic Laplacian. Speaking about this discrete model we do not 
mean just the corresponding difference operator on a lattice or on graphs but we 
mean a discrete analog of the Riemannian structure on some combinatorial object. 
We consider discrete forms as certain cochains. We construct discrete analogs of 
the exterior multiplication operation, the Hodge star operator, of inner product 
(II. 1|) and the operators (|1.2(l . JT3J). 

Our approach bases on the formalism proposed by Dezin PJ. For an account of 
other geometric finite-difference approaches to Hodge theory of harmonic forms see 
references The discrete magnetic Laplacian on graphs had been studied in 

In the spirit of [3], 0, E3 we study self-adjointness of the discrete magnetic 
Laplacian and we proof that the Dirichlet problem for the discrete Poisson type 
equation has a unique solution. 

In this paper we consider just the two dimensional Euclidean case. Although sim- 
ilar constructions can be carried out in the n-dimensional case, the two-dimensional 
discrete model makes it possible to analyze in detail the combinatorial relations and 
the limiting process. One of the formal results is the construction of a nonstandard 
approximation of the generalized solution of the Poisson type equation for the mag- 
netic Laplacian p. 3(1 under the minimal requirements of smoothness of the right 
hand side (is belonged to L 2 ). 
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2. Preliminaries on combinatorial structures 

We use the schema of discretization due to Dczin 0. Let {xk}, {ek}, k £ Z, be 
the sets of basis elements of real linear spaces C°, C 1 . We will regard the linear 
combinations a = ^a k x k , b = ^b x k , a k ,b k £ R, as zero-dimensional and 
one-dimensional chains, respectively. 

It is convenient to introduce shift operators 

rk = k + 1, ak = k — 1 

in the set of indices. We define the one-dimensional complex C as the direct sum 
C° © C 1 of the introduced spaces with the following boundary operator 

dx k =0, de k = x rk —x k , k £ Z. 

The definition of d is linearly extended to arbitrary chains. We call the complex C 
a combinatorial model of the real line. The basis elements x k , e k can be interpreted 
as points and intervals connecting the points (i.e. e k = (x k ,x Tk )) of real line. 

We consider the tensor degree C{n) = (g>™C of the one-dimensional complex C as 
a combinatorial model of R™. The main object of our study will be a discrete model 
of the magnetic Laplacian in the simplest two-dimensional domain. Therefore we 
describe the combinatorial relations that are encountered in the two-dimensional 
case. 

The basis elements of the two-dimensional complex C(2) can be written as follows 

x k ® x s = x k , s , e k ®x s =e\ s , 
e k <3e s = V ktS , x k ®e s = e 2 ks . 
The boundary operator d we define as 

dx k)S = 0, 9e\ s = X Tk)S — X ktS , ^ e fc,s = x k,rs ~ x k , s , 

(2-1) dV k;S = el >a + e 2 TkiS - e\ iTa ~ e\ s . 

Let us introduce an object dual to C(2). Namely, the complex of complex- valued 
functions over C(2). The dual complex K{2) we can consider as the set of complex- 
valued cochains and it has the same structure as C(2), i.e. K(2) = K®K. In other 
words, K(2) is a linear complex space with basis elements {x k ' s , e\' s , e 2 ' s , V k ' s }. 

The pairing (chain-cochain) operation is defined by the rules: 

(2.2) < Xk ,„ a?« >=< ei >a , >=< e 2 fc>s , e%« >=< V k , B , V™ >= 6™, 

where 5^' q s is Kronecker symbol. We call elements of the complex K{2) forms. Then 
the 0-, 1-, 2-forms ip, u> = (u,v), 77 can be written as 

(2.3) <p = J2<Pk, s x k ' s , w = Y,(u k , s e\' s +v k , a e k 2 ' s ), rj = ^ rjk, s V k ' s , 

k,s k,s k,s 

where tp k , s , u k , s , v k _ s , i] ktS £ C for any k, s £ Z. The pairing (|2.2(l is linearly 
extended to forms l|2.3|l . The boundary operation d in C(2) (|2.1|l induces the dual 
operation d c in K(2): 

(2.4) < da, a >=< a, d c a >, 

where a £ C(2), a £ K(2). The coboundary operator d c is a discrete analog of 
the exterior differentiation operator d. We will need the expression for d c over the 
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basis elements of C(2): 

< 4,8! > = frk : s - <Pk,s = Ak<Pk,s, 

(2.5) < e 2 ks , d c ip >= ip ktTS - ip k ,a = A s ^fc, s , 

< V k , s , d C LU >= V TkiS - V k ,a - U k ,TS + U k ,s = &kVk,s - &sU k ,s- 

We define the multiplication U in K(2) by the rules: 

(2.6) x k ' a U e k ' s = e k { s U x Tk ' s = e k ' s , x k > s U e k 2 ' a = e k ' s U x k > TS = e k '% 

x k,S u F fc, S = u ^Tfe.TS = e fe,8 y e Tk,S = y k , S ^ 

supposing the product to be zero in all other cases. To forms l|2.3|) the U-multi- 
plication can be extended linearly. For arbitrary forms a, /3 £ K{2) we have (see 
|3 p. 147]) the relation 

(2.7) d c (aU(3)=d c aUP+{-l) r aUd c (3, 

where r is the dimension of the form a. So the U-multiplication is an analog of the 
exterior multiplication A for differential forms. 

Let e k,s be an arbitrary basis elements of K(2). We introduce the "star" operator 
setting 

(2.8) e k ' s U *e k < s = V k ' s . 
Using l|2.6[) . we have 

The operator * is extended to arbitrary forms by linearity 
Let now 

(2.9) V = Y,Vk,s, k = l,2,...,N, s = l,2,...,M 

k,a 

is some fixed "domain" , namely, a set of 2-dimensional basis elements of C(2). Then 
the relation 

(2.10) (a, (3) v =<V, a U */3 > 

gives a correct dehnition of the inner product for forms of the same degree (cf. 
I|l.l|)). For forms of different degree the product l|2.10|l is equal to zero. Using 
(E3-<E3, we obtain 

(2.11) (a, P) v = J2 a k,s0k~s, 

fe . s 

where a k)Sl (3 kyS are components of the forms a, (3 £ K{2). 

We agree that in what follows, unless the limits of summation are specified, the 
subscripts k, s always run over the set of values indicated in l|2.9|l . 

Taking into account 1)2.7(1 . (|2.1U|) . we can written for a (p — l)-form a £ K(2) 
and p-form (3 £ K(2) the relation 

(2.12) (d c a, p) v =< dV, aU*^> +(a, 5 C /3) V , 
where 

(2.13) 5 c f3= (-l) p *- x d c * /3. 
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Here * _1 is the operation inverse to *, i.e. = 1. If the form a U */3 vanishes 

on the boundary dV, then Equation l|2.13[) defines the formally adjoint operator to 
d°. Let ui = (u,v) be an 1-form (|2.3[) . Then we have 

(2.14) S c lu = J2(- A kU*k, s - A s v k><TS )x k ' 3 . 

k ■ s 

We call the operator 8 C a discrete analog of the codifferential 8. 

Therefore a discrete analog of the Laplace operator can be defined as follows 

-A c = S c d c + d c 6 c . 

If ip is a 0-form, then — A c ip = 5 c d c ip and we obtain at the point Xk, s the difference 
expression 

< Xk,s, ~A C (f >= 4tpk,s ~ Vrk,s - fk.rs - <£crk,s - <£k,as- 

It should be noted that the definition of the inner product l|2.11[) turns the linear 
space of forms over V into finite-dimensional Hilbert spaces H°, H 1 , H 2 with bases 
{x k ' s }, {e k ' s ,e k ' s }, {V k ' s }, k = 1,2, ...,N, s = 1,2,..., M, respectively. Thus we 
can regard the operators d c , S c , A c over V as follow 

d c :H P ^ H p+ \ S c :H P ^ H p -\ A c : H p -> H p , 

where p = 0, 1, 2. It is convenient to suppose that H^ 1 = H 3 = 0. 

3. Discrete model of the magnetic Laplacian 
Let a real- valued 1-form 

A\ ,A\ £ R, be a discrete analog of the magnetic potential. We define the discrete 
analog of the deformed differential (|1.2|) as follows 

(3.1) d c A :H°^H\ <p^d c ip + itpUA. 
Taking into account l|2.5|) . I|2.6[) . we have 

(3.2) d c A (f = {(Akfk.s + itpk,sA\ :S )e k { s + (A s <p ktS + iip k , s A 2 k ^)e k ' s ) . 

k,s 

As in the continual case (see HI. 40 ). wc can identify the discrete magnetic potential 
A with the multiplication operator 

(3.3) A:H°^H\ tp^ipUA. 
Then we have 

A< P = ^2,{Vk,sA\ s e k i S + ¥k, s A 2 k A S )- 

Proposition 3.1. The formally adjoint operator A* : H 1 — ► H acts on an 
arbitrary 1-form oj = (u, v) as follows 

(3.4) A*oj = J2( A lsUk,s + At !S v k , s )x k ' s . 
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Proof. Since the 1-form A £ H 1 is real-valued by assumption, we have 
(Atp, oj) v = (<p\jA, uj) v =< V, (ip U A) U *uJ > 

k , S 

= ^ s ( j4 M Mfc . s + A k,s V k,s) = if, A*Uj) V . 

k , S 

□ 

Let us suppose that components ak, s of an arbitrary r-form a € H r , r = 0,1,2, 
satisfy the following "boundary conditions": 

(3.5) Qo,s = a T N, s = 0, akfi = au, T M = 
for all k = 1, 2, N, s = 1,2, M. 

Proposition 3.2. Let components of <p € H°, uj <G H 1 satisfy Conditions fff. 5}) . 
Then 

(d c A ip, lu) v = (<p, 5 c A u) v , 

where 

(3.6) 5 c a uj = 8 c uj - iA*uj. 

Proof. Note that Conditions i|3.5|) imply the relation < dV, (pLi*uj >= |3 p. 161]. 
Then from (|2.12[1 we have 

(d c (p, uj) v = (tp, S c uj) v . 

Hence 

(d A tp, ui)v — (d c (p + itp U A, ui)v = (d c ip, us)v + i(<p U A, u)v 
= (ip, S c uj) v + i(ip, A*uj)v = (<P, (S c - iA*)uo) v . 

□ 

Thus the operator 8 C A : H 1 — > H° is the formally adjoint operator to the operator 
d c A . Using (|2.14() . we can rewrite 1|3.6|) in a "pointwise" form: 

< Xk, 3 , 5 c a uj >= -A k u akiS - A s v k:CFS - i(Al tS Uk, s + A\ s v k , s )- 

We have 

ip U 8 c uj = ^2ip k:S (-A k u ak:S - A s v k:Crs )x k > s 

k . s 

= Y J { l P l jk,sUak,s ~ tPk,aUk,s ~ ¥k, S Vk,s + l Pk,asVk,as)x k ' S 
k . s 

+ y^(yfc, 3 ^o-fc,s " <P<,k,aV><rk,8 + V>k, s Vk,as ~ ( Pk,a S VkMs)x k ' S 
k.s 

= 5 C (tf (JLJ)+Y^ {{^kfak^Ucrk.s + {A s <p k Ms)vk,*s)x k,S - 
k . s 

It follows that 

S C ((fi U Uj) = (p U S C U> - ^ ((A k (pak, S )Uak,s + ( A s (p k ,crs ) V k ,as ) X k ' S . 
k . s 
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From this we immediately obtain the following discrete Leibniz rule for S A : 

5 A (p U lo) = (5 C - iA*)(tp Uw) = 5 c (tp U u) - iA*(<p U u) 

= ip(J 8 c lo - pU iA*us - 2J ((Ajfe¥>o-fc,a)«<7jfe,« + (A s ¥>fc, CTS )ufc l0 - s )a; fe ' s 

= (/> U (5^a> - ^ ((A fe ^ CT / IiS )ii <7 fe lS + (A s p k ^ s )v kyrJS )x k ' s 

k . s 

(cf. Sect. 2], where the corresponding Leibniz rule is given in the continual case). 
Let us define the discrete magnetic Laplacian as 

-A c A = S c A d c A :H a ->H a . 

Note that we assume that Conditions (|3.5|) are satisfied for any form p e H°. This 
gives us the necessary extension of p beyond H to consider the operator — as 
above. 

Using IpTTJ) . we have 

-A c A tp = 5 L A (d c <p + ipL)A) 

= (5 C - iA*)d c p + (5 C - iA*)(ip U A) 
= -A c p - iA*d c tp + i8 c (p U A) + A*(p U A) 
(3.7) = -A c p-iA*d c p + i8 c Ap + A*Ap. 

Proposition 3.3. The operator —A C A is self-adjoint, i. e. 

(6 A d c A p, ip) v = (p, 6 A d c A ip)v- 

Proof. It is known (see p. 163]) that under Conditions (|3.5|l the discrete Lapla- 
cian — A c = S c d c : H° — > H° is self-adjoint. Using Propositions 13 . II 13 . 21 we obtain 

(S c A d c A ip, ip) v = (5 c d c p, tp) v - (iA*d c p, ip) v + (iS c Ap, tp) v + (A'Aip, rp) v 
= (p, 5 c d c ip) v + (d c p, iAip) v - {Ap, id c ip) v + (Ap, A%p) v 
= {p, -A c ip) v + (p, i5 c Aip) v - {p, iA*d c i') v + (p, A* Aip) v 
= (p, {-A c + i5 c A-iA*d c + A*A)iP) v = (ip, -A c A ip) v . 

a 

Using 1)3. we can write 



(3.8) 



(d A <p, d c A i/j) v = [d c p,d c %p) v + (d c p, iAip) v 
+ (iAp, d c tp) v + {iAp 1 iAip) v . 
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Taking into account Q2.5(l and Q2.11|l . we have 

(d c tp, A*p) v = J2 (Ak<PkA<h^ A k, s ) + A s<Pk, s (*h7sAls)) 

k,s 

= ^2<Pk,s(- A k (^~ s Al k s ) - A s (V^74 iCTS )) 

k . s 

+ X] (^TJV,*(^JV,a^3v,*) ~ ^(^0^0,5)) 

+ X! {Vk,rM{^k,M A l,M) - fkj{^k,oA 2 kfl )) 
k 

+ ( ( Pk,rM('4'k,M A l,M"> ~ PmW'M^q)) + (y, <5 C ^-0)v- 

It follows that 

S 

s 
k 

~ X] [^,1(^,1 ~ Ipkfi + ii>k,0 A kfi)} 
k 

+ (ip,- A c ij.>) v + {ip,iS c Ail;) v - (ip,iA*(d c <ip)) v + {<p,A*Aip) v . 
Consequently, if the forms <p,tfj € H° satisfy Conditions (|3.5[) . then we obtain 

(3.9) (d A ip, d c A ip) v = - ^2 <Pi,8ipi,8 - ^mV'M + (f, ^ A d c A ip) v . 

s k 

Theorem 3.4. For any form f £ H° a solution of the equation 

(3.10) -A° A <p = f 
exists and is unique. 

Proof. By virtue of the sclf-adjointncss of the operator — it is enough to prove 
the uniqueness of the solution. Assume that (p = ip in Equation 1)3.9(1 . Then we can 
write 

(3.11) (d c A <p, d A tp) v + J2 l^-l 2 + E l^il 2 = ~ A a^)v- 

s k 

Using (E3J), l|2Tr)l . we get 

(<fV, d»y = X (|A fc(/ 3 M | 2 + |A s(/ 3 M | 2 ), 
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(d c ip, iAip) v + (iA<p, d c ip) v = ^E [Ak,s(<Pk,i,(Ak<Pk,s) - (Akfk^wX)] 

k , S 

+ i E [ A k,<,(<Pk,s( A »<Pk,s) - (A a tp kiS )ipk^)], 

k . s 

(iAtp, iA<p) v = E (« s ) 2 |^, s | 2 + (4,J 2 |^ >S | 2 )- 

k . s 

It is easy to check that 

tpk <s (Ak<Pk,s) - {Ak<Pk,s)W^ = 2i(lm(ipjb iS )Re(A fc (p feia ) - Re(^fc, a )Im(A fc <p fcia )). 
Substituting the last relations into i|3.8|) we obtain 

{d^p, d c A p) v = E [l A ^Ml 2 + |A S ^ M | 2 + (4, s ) 2 |^,.| 2 + (A% s ) 2 \<p k , s \ 2 

k,s 

+2Al s Re(<p ktS ) hn(A k ip ktB ) - 2A x k s Im(tp k>s ) Re(A k <p ktS ) 
+2A 2 ks Re(tp k , 8 ) Im(A s v? fc , s ) - 2A 2 ks Im(<p k , s )Re(A s (p ktS ) 

= 5Z [( Re ( A kfk, s ) - Al s Im{ip Ks )y + (lm(A fe 95 fe . s ) + ^ s Re(v3 fc , s )) 2 

+ (*Re(A s ^fc, s ) - Al iS hn(ip k ,»)) + (*Im(A sVfc , s ) + Re(ip k>s f) . 

Now let we take / = in Equation H3.10|l . Then comparing the last equation 
and l|3.11|) . we obtain 

E [(Rc(A^ fc , s )-^Llm(^fe. s )) 2 + (lm(A fc (^ M )+4 s Rc(^ M )) 2 
+ (Re(A s¥ > M ) - A\ s Im(93 fejS )) + (lm(A s p M ) + vl^ s Refto*,,)) 

+E!^mI 2 + EI^mI 2 -o. 

s fc 

It follows that c/3/t.s = for any k, s. Hence vp = 0. □ 

This immediately implies the following statement: 
Corollary 3.5. TTie operator —A C A is positiv. 

4. Approximation and limiting process 

In this section we consider the relationship between the combinatorial objects 
that we have described above and the corresponding continual objects. We will con- 
struct some nonstandard approximation of the generalized solution of the Poisson 
type equation 

(4.1) -A A p = f, 

where / 6 L 2 (fl). We will realize the scheme similar to that given in Ch.3, 
Sec.3]. 
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Let the domain O G R 2 be a rectangle with vertices (a%, bi), (02, 61), 
(ffli, fe), (02, 62), where < ai < 02, < fei < 62- We introduce a scale h setting 
h = N~ 1 (a,2 — 0,1) = A/ _1 (&2 — bx). Divide f2 by the following straight lines 

x = 01 + fc/i, y = &i + s/i, fc = 0, 1, N, s = 0, 1, M. 

We denote by xt. s the point of intersection of these lines. We denote by V& )S an open 
square bounded by the lines: x = a\ + kh, y = b\ + sh, x = a\ + rkh, y = b\ + rsh. 
Let e\ s and e 2 s be the horizontal and vertical sides of V^ s , i. e. e\ s = (xk^:X T k,s): 
e t s = ( X M> lfc,Ts)' I n this way we identify the rectangle ft with the combinatorial 
domain V lfO|) . 

Let us now compare every discrete form ip G H with the step function assuming 
that 

ip h (x, y ) = tp k , s , fox x,y eV k ,s- 
In the case of the 1-form u> = (u, v) E H 1 we have the pair of step func- 
tions u h (x,y) = Uk tS , v h (x,y) = Vk, s and we can write u> h (,x,y) = u h (x,y)dx + 
v h (x,y)dy. Recall that <Pk,s,Uk, s ,Vk.s & C for any k,s. Then p h ,uo h are complex- 
valued. 

It is easy to check that 

(4-2) = fr-IMIffo, II^Hz^AMn) = h\\u\\m- 

Define difference operators acting on the step functions as follows 
A h x <p h {x, y) = h- 1 ( 9 h (x + h,y)- <p\x, y)) , 
A h y p\x,y) = h- x {ip h {x,y + h) - V \x,y)). 

Replacing the partial derivatives Jj, appearing in d, S by the difference op- 
erators Aj, Ay, we can introduce the difference operators d h , S h . The difference 
equation 

(4.3) d h ip h = uj h 

is equivalent to the following family of equations 

A k ip ki s = hu k , s , A s ip kiS = hv k , s , 

where k = 0, 1, N, s = 0, 1, M. Hence Equation (|4.3|) can be rewritten as the 
following discrete equation 

d c p> = huj, 

where ip, to are discrete forms (see (|2.3|) ) with the components (fk. s and Ufc, s , i>fc lS , 
respectively. Similarly, we associate the difference equation S h oj h = p h and the 
discrete equation 5 c uo = hip. 

We can also introduce the following difference operator 

(4.4) d h Ah =d h + iA h , 

where A h — A lh dx + A 2h dy and A lh ,A 2h are real- valued step functions defined as 
above. On the other hand, we can consider the step 1-form A h as the multiplication 
operator A h : L 2 (Q) — > L 2 A 1 (fJ) acting as follows 

A h p h = p h A lh dx + p h A 2h dy. 

Then the formally adjoint operator to A h (cf. I|3.4[) ) acts on a step 1-form oj h = 
(u h , v h ) as 

(A h )*uj h = u h A lh +v h A 2h . 
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Thus wc define the difference magnetic Laplacian (cf. I|3.7jl ) by the formula 

(4.5) -A h Ah = 5 h Ah d h Ah = 5 h d h - i(A h )*d h + i5 h A h + (A h )*A h . 

Now we consider the discretization procedure (see for details p. 170]). Let 
f(x,y) be a complex-valued function defined over Q (or over V = Vk,s) and let 
/ £ L 2 (Q). Associate / with the step function f h setting 

(4.6) f h (x,y) = h- 2 J f(£,n)d£dti, for x,y€V h , s . 

v k , s 

Moreover, the value of f h can be assigned to the point Xk.s- As above, we can write 
f h (x,y) = fk, s for x, y £ T4. s . Thus we obtain the discrete 0-form / = fk, s x h ' s 
which is associating with / £ L 2 (£l). Similarly, if / is an 1-form, / £ L 2 A 1 (il), 
then we associate each component of / with the step function (j4.f)|l and we assign 
the value of it to one of the intervals e\ or e^ s . 
Let us introduce the norm 

|| <p h f w{Q) = J (\A h x ^\ 2 + \A h y\ 2 )dxdy. 

n 

It is not difficult to verify that 

N M 

(4-7) II ^ ||^ r ^^(|A^ M | 2 + |A^ M | 2 ). 

fe=0 s=0 

Hence we can write 

II V h \\w{n) = \\ V \\w{v) ■ 

Theorem 4.1. Let the step function f h be the discretization of f £ L 2 (Q). Then 
the following Dirichlet problem 

(4.8) ~A h Ah ip h = f\ 

(4.9) </|an = 
has a unique solution and the inequality 

(4.10) || ip h || W r(n)<ci||Rfi/|| iB( n)+Pa||lDa/||£2(n) 
is valid for the solution ip h . 

Proof. Using l|4.5|l . Equation l|4.8|l can be rewritten as 

(4.11) 5 h d h ip h - i{A h )*d h Lp h + i8 h A h tp h + (A h )*A h <f h = f h . 

By definition the step function ip h and the step form A h on fl associated with the 
discrete forms ip and A on V. As above, if we replace the difference operator d h , 5 h 
by the discrete operator d c ,S c , then Equation (|4.11|) transforms into the following 
equation 

(4.12) S c d c <p - ihA*d c ip + ih6 c A<p + h 2 A*A<p = h 2 f, 

where / is a 0-form defined by the step function f h (see H4.60 ). Note that, if the 
step function ip h satisfies Condition (|4.9|) . then the corresponding discrete form 
satisfies Condition l|3.5|) . Thus the unique solvability of Q4.8[l . (|4.9|) immediately 
follows from Theorem 13.41 
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Let now represent Equation (|4.8|l as follows 

~A h Ah Re ip h = Re f\ -A\ h Im tp h = Im f h . 

In a similar way we can split Equation (|4.12(l . 
Since 

(d c a, iAa)v + {iAa, d c a)y = 
for any real- valued discrete form a £ H° it follows that from i|3.7|) , i|3.11|l we obtain 

II^Re^H^ + ^(Rc<p hs ) 2 + ^2(Re(p kA ) 2 + h 2 \\ARc<p\\ 2 H1 = h 2 {Rc<p, Rcf) v . 

s k 

It immediately follows that 
II R e vllv7(y) 

<h 2 (Rcip, Re/) v = /i 2 ^Rc^ M Rc/ M . 

fc,S 

Hence, using (|4.2(1 . we have 

(4.13) || Re/||^ (n) < || Re/|| L2(0) || Ref h \\ L 2 (a) . 
It is easy to check the following estimates 

||Re^||i2 (n) < || Re ^H^n), II Ref h \\ LH n) < Cl || Rc /|| L2(n) 
(see for details jS] Ch.3, Theorem 5]). Combining the last with (|4.13|) . we obtain 

(4.14) \\Rep h \\ wm < Cl ||Re/|| L2(n) . 
Similarly we obtain the estimate 

(4-15) || Im<p h \\ w(Q) < c 2 || Im/|| L 2 (0) . 

Finally since 

||^' l ||w'(o) < II Rctfi h \\w(n) + II Im-P h \\w(n) 
Estimates lf4"TIj) . igUSjl imply QETty . □ 

Let us now consider the limiting process. As in [3| Ch.3], by the step function 
ip h we construct the smooth (i.e. of the class C 1 ) function J h Lp h . It is convenient 
to take J h if h in the form 

x+h y+h 

J h ^ h (x,y) = h- 2 f [ ip h ^,r,)d^dr,. 



This is the well known Steklov function with the averaging radius equal to the 



parameter h (the scale of the net). We can also define the 1-form J h A h £ AhJQ) 



as 

J h A h (x,y) = J h A lh (x,y)dx + J h A 2h (x,y)dy. 

Denote by W 1 ^) the Sobolev space of complex-valued function which satisfy 
the homogeneous Dirichlet condition. 

Let us consider some sequence {h n } such that h n — > as n — > +oo. For conve- 
nient further we will write h instead h n . 

Let A G Ajj^ri) in the operator — A^. We have the statement. 
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Theorem 4.2. Let a step function ip h be the solution of the Dirichlet problem {4-8( 1, 
14-9(1 for the given element f E L 2 (fl). Then the sequence {<p h } strongly converges 
in L 2 (fi) to the element ip £ W 1 ^) as h — > ; where ip is the generalized solution 
of the corresponding Dirichlet problem for Equation \4-l\l - At the same time the 
sequence {J h (p h } converges to ip in the metric W' 1 (J7). 

Proof. Based on Theorem l4.ll the proof is similar to that of Theorem 5 ch.3] . □ 
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